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We study the Kondo lattice model using a class of canonical transformations that allow us to faithfully rep- 
resent the model entirely in terms of fermions without constraints. The transformations generate interacting 
theories that we study using mean field theory. Of particular interest is a new manifestly 0(3)-symmetric rep- 
resentation in terms of Majorana fermions at half-filling on bipartite lattices. This representation suggests a 
natural 0(3)-symmetric trial state that is investigated and characterized as a gapped spin liquid. 



I. INTRODUCTION 

One of the standard models in strongly correlated con- 
densed matter physics is the Kondo lattice model, for reviews 
see e.g. Refs. 1^. The basic physics that this model aims to 
describe is the simultaneous existence of and interaction be- 
tween conduction c-electrons and localized /-spins. It is chal- 
lenging to describe both subsystems on equal footing because 
the algebras of the creation and annihilation operators of the 
c-electrons and the /-spin operators are quite different. The 
aim of this work is to introduce and study faithful fermionic 
representations of this model. 

There exists many different approximation schemes 
that have been used to study this model, for example 
Gutzwiller projection,^-^ slave particle methods,^ large-N 
approximations,***^ canonical transformations,"' mean field 
theories," '- and dynamical mean field theory.'^ For discus- 
sions of these methods and more references we refer to the re- 
views just cited. Numerous numerical studies have also been 
performed. These are particularly successful in one dimen- 
sion (ID),^ but also the two-dimensional (2D) Kondo lattice 
have been studied using quantum Monte Carlo at half-filling 
where the sign problem is absent. 

All approximate treatments have their shortcomings. The 
validity of large-N approximations are questionable when one 
considers the physical case of = 2 for example. Slave- 
particle representations are popular because they keep the spin 
and charge symmetries clearly visible. One drawback of this 
approach is that the exact constraint that should be imposed is 
typically only imposed on average and not exactly. This can 
be improved at the expense of introducing gauge fields.'^ In 
this work we use exact canonical transformations and hence 
this drawback will not affect us. 

Another physically appealing approximation scheme in- 
volves a simple mean field decoupling, see e.g. Ref. 16. This 
is found to be in qualitative agreement when compared to 
more involved numerical methods.'^ One problem with this 
approximation is that it does not reproduce the correct size 
of the Hilbert space of the Kondo lattice. Because the mean 
field Hamiltonian has two bands for each spin component, the 
size of the Hilbert space is 2^ per unit cell. This is appro- 
priate for the periodic Anderson model when the onsite inter- 
action strength U is small. In the Kondo lattice there are 2"^ 
states per unit cell however. This problem will not affect us 
in this work since we will construct faithful representations 
using three fermions per unit cell. 



As we will see there is considerable freedom in the allowed 
canonical transformations. Of particular interest is a new 
manifestly 0(3)-symmetric representation on bipartite lattices 
at half-filling. This representation is most easily described 
in terms of Majorana fermions. One way to motivate it is to 
start from the known representation of a spin-^ in terms of 
three species of Majorana fermions: Sj- = —ie^^^fijUk/^. 
This representation has a long history that goes back to the 
fifties, see e.g. the references in Ref. 18. Similar represen- 
tations have become increasingly popular in recent years af- 
ter the introduction of the Kitaev model. This representation 
was also heavily used in earlier works focusing on non-Fermi 
liquid behavior in modified Kondo impurity problems,^"'^' as 
well as lattice systems,^" and odd-frequency pairing in Kondo 
lattice models. The fact that this representation can be used 
to study Heisenberg models and generate spin liquid states 
has also been known for some time.'^'^^-^"* A variant of such 
a spin liquid state on the triangular lattice was suggested very 
recently.'^ 

It is also well appreciated that the representation of a spin- 
I in terms of three Majorana fermions is redundant. In fact 
the operator 70 = 2i/xi/Lt2M3 commutes with the spin operator 
and can be viewed as an independent Majorana fermion. The 
0(3) representation can be obtained by writing the c-electron 
creation and annihilation operators in terms of three other in- 
dependent Majorana fermions and this composite Majorana 
fermion 70 . The Heisenberg exchange interaction between the 
/-spin and the spin of the c-electrons then takes on a simple 
form with manifest 0(3)-symmetry. To the best of our knowl- 
edge this procedure has not been used previously to describe 
the standard Kondo lattice model. 

This construction can be straightforwardly generalized to 
the Kondo lattice. But the local 0(3) symmetry is a result 
of the underlying SO(4) symmetry of the Hubbard model at 
half filling on bipartite lattices,-^ and its generators are com- 
binations of spin and pseudospin generators (see Sec. III). It 
is therefore only possible to generalize the local 0(3) symme- 
try to a globally 0(3)-symmetric Hamiltonian on bipartite lat- 
tices. On bipartite lattices it is then easy to write down a sim- 
ple translationally invariant 0(3)-symmetric mean field state. 
This state is found to be favorable when the kinetic term is 
small or moderately large compared to the exchange interac- 
tion. The resulting state is a gapped spin liquid with nonzero 
triplet pairing amplitude. 

The paper is organized as follows: For completeness we 
end this introduction by writing down the standard form of 
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the Kondo lattice Hamiltonian explicitly. In Sec. II we dis- 
cuss canonical transformations for the 1-site Anderson model 
and derive different representations for this model. In partic- 
ular we consider the limit of this model that reproduces the 1- 
site Kondo impurity problem. These local transformations can 
easily be generalized to the lattice since they involve canonical 
fermions. In Sec. Ill we provide an elementary derivation of 
the 0(3)-symmetric representation of Sec. II using Majorana 
fermions. We also consider the hopping of c-electrons. On bi- 
partite lattices the hopping term can be chosen in a particularly 
symmetric form with a manifest global 0(3)-symmetry. The 
resulting Hamiltonian is written down explicitly in Eq. (21) of 
Sec. IV where we study this interacting fermion model using 
variational mean field theory at half-filling. The variational 
solutions are worked out in ID for simplicity. Of particular 
interest is an 0(3)-symmetric mean field state that is found to 
be favorable for some parameter values. This spin liquid state 
is characterized further in Sec. V where it is shown to have ro- 
tationally invariant spin-spin correlations as well as a nonzero 
triplet pairing amplitude. Our conclusions and an outlook are 
to be found in Sec. VI, and some mathematical details in an 
Appendix. 



A. The Kondo lattice model 

In this paper we study the Kondo lattice model that can be 
described by the Hamiltonian- 

Hklm = -tJ2Yl 4,a(n)cc,^(r'j) +h.c. 

-M^K(r.)-l] + J^5,(r,)-5/(n). (1) 

i i 

The first line describes the hopping of c-electrons on the lat- 
tice, and (T = i is a spin label. For simplicity we consider 
only hopping between nearest neighbor lattice sites and rj, 
as indicated by the notation j) ■ The first term on the sec- 
ond line can change the average number of c-electrons away 
from half-filling by introducing a non-zero chemical poten- 
tial fi. ncivi) = J2a ^ta{i^i)Cc,a{i^i) the number opera- 
tor of the c-electrons. The second term on the second line 
describes the interaction between the spins of the c-electrons 
and the localized /-electron spins. On each site of the lat- 
tice, Sf{ri) is a local spin-| object (we use units such that 
h — 1 throughout the paper) that satisfies the usual spin- 
i algebra [S\ S^] = ie'^'^S'' and = 3/4.2^ The spin 
operator of the c-electrons can be represented as Sc{ri) — 
\ Y.a,a' 4,<T(^i)'^<T,<T'Cc,o-'(n), with t'' (a = 1 ... 3) the stan- 
dard Pauli matrices. 



II. CANONICAL TRANSFORMATIONS FOR THE 1-SITE 
ANDERSON MODEL 

The spin-spin interaction term in (1) can be derived from 
the periodic Anderson model in the limit that the occupation 



of the /-electrons is not fluctuating. For this derivation it is 
sufficient to consider a 1-site Anderson model on each site of 
the lattice. This is the model that we consider in this section 
and we extend the formalism in Ref. 10 to allow for more gen- 
eral canonical transformations. Because we are dealing with 
only one site we will suppress the site index in this section. 



A. 1-site Anderson model 

We consider first the symmetric 1-site Anderson Hamilto- 
nian at half filling 

^fiA = -WY. {cl^Cf^, + c]^^c,J + U{nf - If, (2) 

where riy — J^a^f-r^ ~ '^a'^f a'^s a- The Hilbert space 
of this model consists of 16 states that are enumerated in 
Eq. (Al). Following Ref. 10 we then introduce another basis 
using electron e\ and hole h\ creation operators, the corre- 
sponding states are enumerated in Eq. (A7). 

We define the states in the new basis according to the fol- 
lowing rules: 1) The singlet ground state |0)s of Eq. (2) is 
equal to the vacuum state \l)eh that has no quasi-particles 
(i.e. no electrons or holes). 2) The lowest energy eigenstates 
with the appropriate charge and spin (measured with respect 
to the singlet ground state) are mapped to the states with one 
quasi-particle, i.e. /i||l)eh, e\\l)eh, h\.\l)eh, and e\\l)eh- 3) 
The assignment of the other states is fixed by the charge and 
spin quantum numbers, except for 4) the two remaining sin- 
glets that can be assigned by considering the parity transfor- 
mation that exchanges the /- and c-electrons. This uniquely 
defines the quasi-particle operators up to a gauge choice. We 
fix this arbitrariness by demanding that the quasi-particles go 
smoothly to the bonding and antibonding eigenstates in the 
non-interacting limit [/ = 0. The transformation matrix that 
implements this transformation is given explicitly in Eq. (A8). 

This transformation is similar to the one in Ref. 10, ex- 
cept that it provides an adiabatic connection between the low- 
energy states in the non-interacting and the interacting sys- 
tems. Since the quasi-particle excitations also have the same 
spin and charge quantum numbers as the non-interacting par- 
ticles this transformation is an explicit demonstration of the 
Landau quasi-particle concept. 



B. Kondo limit ?7 ^ oo, fermionic representations 

The Kondo limit can be obtained by taking J7, — > oo 
while keeping J — AW'^/U finite. In this limit the Hamilto- 
nian (2) becomes 

HiA Hj + U{nf ~ If = JSc ■ Sf + U{nf - if. (3) 

The last term is effectively a constraint that enforces that the 
/-level is singly occupied. In the Kondo limit it is therefore 
natural to describe the remaining low-energy spin degree of 
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freedom of the /-electrons in terms of a localized spin- ^ op- 
erator Sf. Doing this on every site of the lattice we are left 
with the spin-spin interaction term in Eq. (1). 

Alternatively this limit can be studied in terms of the elec- 
tron and hole operators of the last subsection.'" This represen- 
tation has the advantage of having nice transformation prop- 
erties of charge and spin. A disadvantage is that states with 
more than one quasi-particle present (outside the triplet sec- 
tor) at a site will have weight in the high-energy sector. This 
weight can nevertheless be kept small in the Kondo lattice by 
being close to the atomic limit t J.'" 

We will instead use a less symmetric representation of 
charge and spin with the basis defined in (All). The advan- 
tage of this is that it will allow us to entirely get rid of the 
high-energy sector by taking one of the operators that span 
the local Hilbert space (i.e. c\) to encode the high-energy ex- 
citations. This operator is then straightforwardly eliminated 
from the periodic Anderson model to generate a low-energy 
theory that is the Kondo lattice model. 

There is considerable freedom in how to assign the 8 
states in the low energy sector to the low-energy states in the 
electron-hole basis. In this work we will only consider trans- 
formations such that the operators cjj (a = 1 ... 3) are odd 
with respect to fermion parity. This has the advantage that 
these operators will be fermionic also in between sites in the 
Kondo lattice. This implies that the low-energy theory is rep- 
resented by a theory of interacting fermions on a lattice that 
can be studied with standard many-body methods. The main 
disadvantage of this representation is that it is not possible 
to do this without mixing the spin and charge content of the 
theory. It is also possible to take one or two of the cj,'s to 
be hard-core bosons, but we will not consider these types of 
transformations in this work. 

With this constraint on the allowed transformations we have 
that; 1) The vacuum with energy Eg = —3 J/4 is the ground 
state singlet \0)s- 2) The charged states with one 

electron or one hole and energy Ec = maps to the states 
with one quasi-particle (|2)x, |3)k, and \4:)x) or three quasi- 
particles (|5)a')- 3) The three triplets with energy Et — J/4 
maps to the three states with two quasi-particles i\6)K, \'^)k, 
and \8)k)- The Hamiltonian (3) in this basis is therefore al- 
ways of the form 

Hj='^[l-{n,+n2 + n3~2f], (4) 

independently of how the states are actually assigned to the 
three fermions. This Hamiltonian is clearly invariant under 
a U(3) rotation in the vector space spanned by the three op- 
erators {c[,C2,cl), but we will restrict ourselves to the real 
subgroup 0(3) in this paper. 



C. At half -filling 

At half filling it is desirable to use a representation that 
is symmetric between electrons and holes, and also as rota- 
tionally symmetric with respect to spin as possible. We have 



found that the following procedure gives quite symmetric rep- 
resentations: 1) We would like the state with three quasi- 
particles to be a state with zero charge and spin on average. 
A possible choice is \5)k = S5(e| — h!^)\l) eh / V^, with 
S5 = ±1. This couples the charge and spin content of the 
excitations, but this is unavoidable in this approach. Note 
that we have arbitrarily picked a direction of the spin of the 
electron part of the state. 2) The other three odd fermion 
parity states can be assigned by demanding that the oper- 
ators c\ and C3 are related by time-reversal symmetry, we 
choose the convention Tc\T^^ — —c^ and Tc^T^^ — c[. 
A possible assignment of the states is therefore \2)k — 
~ie\ - hl)\l)e,jV2, \3)k = (4 + hl)\l)eh/V2, and 
|4)a' = (e| + h\.)\l) eh / V^- 3) Because c\ and 4 re- 
lated by time-reversal symmetry T| 7) ^ — \7)k. This implies 
that \7)k = -S7(|6)eh + |8)e/,.)/\/2, with 57 = ±1, since this 
is the only time-reversal invariant triplet (see the discussion in 
Sec. A 2). 4) The assignment of the remaining two triplets can 
be parametrized by an angle 0t and a sign sq — ±1 according 
to Se\6)K = cos{(j)t){\6)eh - \8}eh)/V2- sin{(j)t)\7)eh and 
\8)k = sm{(j>t){\6)eh - \8)eh)/V2 + cos(0t)|7>e/>. 

The generic transformation is therefore parametrized by the 
numbers (55 sq S7,(j)t)- Working out the transformation and 
expanding the original c-electron operators in terms of the new 
fermions we find that they are typically of fifth order This ex- 
pansion is most easily performed with the aid of a computer 

Only the combinations (+++, 0), (— + — , 0), and ( h, ipt) 

terminate at third order The (+ + +, 0) transformation gen- 
erates the most symmetric representation: 

J _ ci + 4 , C2 4 
22' 

t C3 + 4 , (ci-4)(c2 + 4)(c3-4) 

Cc,i- 2 + 2 - ^' 

This representation is most easily formulated in terms of Ma- 
jorana fermions. An alternative elementary derivation of this 
representation is provided in Sec. III. The other transforma- 
tions are equivalent up to a rotation, as an example we give 
the expression for the representation generated by ( h, tt): 

^ _ (ei+4)(l-2n2) , (c2-4)(l-2ni) 
"^-t- 2 + 2 ' 

t _ C3 + 4 , (ci + 4)(c2-4)(c3-4) 

Cc,i- 2 + 2 ■ ^' 

D. Away from half-filling 

Away from half-filling there is no reason to try to enforce a 
symmetry between electron and hole excitations. Considering 
the case of hole doping the hole excitations will have lower 
energy than the electron ones. It is therefore natural to choose 
two of the low-energy creation operators (c| and Cg say) to 
create the two hole states. With this choice the chemical po- 
tential term becomes < for hole doping) 

= -^i{nc - 1) = -^J-{n2 - ni - 713 + Snins). (7) 
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To be concrete we define \2)k = ^||0)s and |4)if = ft.||0)s. 

It is also convenient to let the creation operators c| and Cg 
to be related to each other by time-reversal symmetry as in 
the case above. This means that the states \Q)k, |7)_r-, and 
|8)if can be parametrized exactly as in Sec. II C. The electron 
states are defined via |3)x = [cos(</ie)e| + sin(0e)e|] \0)s 

ands5|5)K = [— sin(0e)e| +cos((/)e)e|] |0)s. Expanding the 
c-electron operators in terms of the new fermions only two 
classes of transformations terminate at third order In both 
cases (j)e = (f>t and the sign structure is (+ + +) or ( — \ — ). 
The representations are (0 = (pt and s — 57) 

t _ - "-3) + gc|"3 



V2 

-^[cos{(f>) + ssin(</))(ci - scl)(c3 - s4)], 

23(1 — ni) + SCgTT-i 

71 



- ^[sin((?!)) - scos{(j)){ci - sc|)(c3 - scl)]. (8) 

We leave the investigation of the representations in Eqs. (6) 
and (8) for a later study and will in the following focus on the 
most symmetric representation of Eq. (5). 



III. THE 0(3)-SYMMETRIC REPRESENTATION AT 
HALF-FILLING 



independent fermion one will somehow have to deal with the 
fact that the Hilbert space has been enlarged, see e.g. the dis- 
cussions in Refs. 24 and 25. The model with an independent 
70 was studied in the context of odd-frequency pairing at the 
mean field level in Ref. 22. The effect of the enlarged Hilbert 
space then shows up as an additional term in the mean field 
Hamiltonian, i.e. Hq in their Eq. (3.3). The remaining part 
of their mean field Hamiltonian is similar to ours, the main 
difference being that they do not generate nonlocal hopping 
terms involving the /i's. 

The c-electron spin operators are 



-^(7273 + 7i7o)/2, 
-«(737i + 727o)/2, 
-«(7i72 + 737o)/2. 



(12) 



In terms of the parity operators pa = 2i7a/ia (a = 1 ... 3), 
which each has eigenvalues ±1, the exchange term Hj = 
JSc ■ Sf can be worked out to be 

Hj = '^{pi +P2 +P3) - '^{piP2 +P2P3 +P3P1) 

= ^{^~PlP2P3){Pl+P2+P3)- (13) 



With the identification pa = 2na — 1 we see that (4) and (13) 
are equivalent. Note that the spin-spin exchange term has be- 
come partly quadratic in the fermions in this representation. 
This is not the case if one treats 70 as an independent Majo- 
rana fermion. 



In this section we will consider the representation in (5) 
from another point of view. As discussed in the introduc- 
tion, it is well-known that it is possible to represent a spin- ^ 
operator in terms of three species of Majorana fermions /i^, 
a = 1 ... 3. These Majorana fermions are real /ijj = 
and independent {/Iq, /ih} = Sab- It is then straightforward to 
check that the operators 

Sj = -iAi2Ai3, = -i^J.3Pl, = -ipiiJL2, (9) 

satisfy the angular momentum algebra [S'j, = ie^'^^'S^ for 
a spin-i since S'j — 3/4. Let us now define 

70 = 2i/xi/X2Ai3- (10) 

This is a proper Majorana fermion (with 7q = 1 /2) that com- 
mutes with the spin operator of the /-electrons: [70, Sj] ~ 0. 
Therefore we can represent another set of Dirac fermion oper- 
ators (the c-electrons) in terms of 70 and three other Majorana 
fermions 7a, a = 1 ... 3, as 

7i-*72 -73 -«7o 

This representation can be found in e.g. Ref. 22, but the nov- 
elty here is to use the composite operator in (10) for 70, in- 
stead of an independent Majorana fermion. Both choices sat- 
isfy the correct operator algebra, but if one keeps 70 as an 



A. Pseudospin symmetry at half -filling 

At half-filling the bipartite Hubbard model, and hence the 
bipartite symmetric Anderson model possess another sym- 
metry. This is called pseudospin symmetry and is imple- 
mented by exchanging the roles of electrons and holes in one 
of the spin components. We can implement this by taking 
7o ^ —70- The generators of the pseudospin algebra are then 



Il = 


-i(7273 


- 7i7o)/2, 




-i(737i 


- 7270) /2, 


Il = 


-i(7i72 


- 737o)/2. 



It is straightforward to check that [P.P] — ie^^^I^ and 
[/*, S^] — 0. If we enforce unit occupancy for the /-electrons 
we see that If = Q so that the pseudospin algebra of the /- 
electrons becomes trivial. Combining the spin and pseudospin 
symmetry the system has a global SO(4) symmetry,-*" which 
is very transparent in the Majorana representation.^*^ 

Going away from half-filling only the generator com- 
mutes with the Hamiltonian and we recognize as the gen- 
erator of the U(l) gauge symmetry related to charge conserva- 
tion. This can be seen by including a chemical potential term 
for the c-electrons in the Hamiltonian: 

Hf, = -fJ.{nc - 1) = i(7i72 - 7370)/^ = -2m/c ■ (15) 
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B. Arbitrary Kondo lattice 

Let us now consider another site with the same representa- 
tion: 



7i - «72 ~ -73 - «7o 



V2 ' ' V2 ' 

The hopping term between neighboring sites then becomes 



(16) 



(7370 + 7370 + 7271 + 727i) • (17) 

This representation has the advantage that it generates at 
most quartic fermion terms in the Hamihonian. This impHes 
that the interaction terms can be decoupled using standard 
Hubbard-Stratonovich transformations. 



C. Bipartite Kondo lattice at half -filling - an 0(3)-symmetric 
representation 

On a bipartite lattice it is useful to employ different repre- 
sentations on the two sublattices. We introduce an extra phase 
of 7r/2 on one sublattice such that 



/ 7i ^ ^72 \ ~ _ . / -73 - ^70 \ 



V2 



(18) 



Because this just involves a gauge transformation it leads to 
the same on-site exchange term. Using this the hopping term 
between two neighboring sites becomes 



fermion theory. For concreteness we will only solve the 
mean field equations for the ID lattice with nearest neigh- 
bor hopping, but the generalization to other bipartite lattices 
is straightforward. The Kondo lattice Hamiltonian in the rep- 
resentation of Sec. Ill C is 

3 

HkLM =-t^Yl *7a(T-07a(rj) 

i j 

- M 5Z [*72 [rihi (n) + (n)7o (r,)] 

i 

- AiXl[*^2('^j)7i('^j) + «73(ri)7o(rj)] • (21) 



Hj{ri) and Hj{rj) are the generalizations of Eq. (13) to in- 
clude a lattice index. The two sublattices are distinguished by 
the absence or the presence of a tilde. 



A. 0(3)-symmetric mean field at half -filling 

The simplest mean field Hamiltonian at half-filling (i.e. 
/i = 0) is manifestly 0(3) invariant: 



3 

^^0(3) — ^ Ha , 
a=l 



(22) 



Ha = ~tY^ nairi)ja{rj) + ^ g{rj ~ r,)ifj.a{r,) flair j) 
+ V^i'-/a{r,)fia{rt) + V^i%{rj)lla{rj)- (23) 



H = -t (c' c 



h.c.) = -i}_^z7„7,. (19) 



In the representation in terms of 7^ and the /i^ (a = 1 ... 3) 
*7o7o is a non-linear operator This implies that we have made 
a non-linear transformation that preserves a fraction of 3/4 of 
the linearity of the hopping term. Explicitly the remaining 
non-linear part of the hopping term can be written as 



«7o7o = -4(iAiiAi)(«Ai2A2)(jM3A3)- 



(20) 



Clearly the Hamiltonian has a global 0(3) symmetry that is 
obtained by the rotation of all local vectors of (71 , 72, 73) and 
(Mi : M2 : M3 ) in the same way. This corresponds to a rotation of 
spin and pseudospin with the same angle, i.e., the three sym- 
metry generators are S"^ + 5j- + for j = 1 ... 3. Note also 
the peculiar feature that the c-electron charge and spin opera- 
tors are no longer quadratic in fermions in this representation. 



IV. MEAN FIELD STUDY 

In this section we will use a mean field Hamiltonian, which 
is quadratic in the fermions, to approximate the interacting 



This form can be motivated from a mean field decoupling 
of (21) using (13) and (20). We now take the ground state 
of the mean field Hamiltonian in (22) as a trial state to ap- 
proximate the ground state of the full interacting theory of 
Eq. (21). Rather than fixing the variational parameters by 
the usual Hartree-Fock decoupling procedure directly, we will 
keep them arbitrary for the time being, since it is in general 
possible that different mean field Hamiltonians give the same 
trial state. 

It is straightforward to diagonalize this problem by going to 
Fourier space using 

7a(r,) = ^^^'[e''=-^7jfc) + e-*'^-7l(fc)], (24) 

and similarly for the other operators. The prime indicates 
that one should only include one of the states for each pair 
of k and — fc in the sum, see e.g. the discussion in Refs. 21 
and 22. This is a consequence of 7q(— fc) — lX{k). The 
physics is independent of the choice of k or — fc. Note 
also that the Brillouin zone corresponds to a lattice with two 
sites per unit cell, hence the N/2, where N denotes the to- 
tal number of lattice sites. Introducing the spinors (fc) = 
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[^a(fc),7a(fc), Aa(fc),7a(fc)]'^ One has 

^^'vI/t(fc)H(fc)vI/Jfc), 



0(3) 



(25) 



a— 1 fc 



where the Hamiltonian matrix 1-L{k) is 



H{k) = 



V 





iV 
-zg*{k) 


ik-S 



-iV ig{k) \ 

-itaijt 

-if 

ita'ik) iV / 



(26) 



Here a{k) = J2j ' and g{k) = gi.^M'^ \ where 
5j are the vectors that connects one lattice site (without a 
tilde) to its nearest neighbors. The mean field Hamiltonian 
is straightforwardly diagonalized in any dimension, but for 
simplicity we only perform the mean field analysis in ID 
in the following. Setting the nearest neighbor distance to 1 
sums can then be converted to integrals with the replacement 
^E' ^ S^'^ t- We now introduce g± = [5(l)±g(-l)]/2, 
and the mean field solution involves solving for g±, V, and V. 
For each value of k the spectrum of'H{k) is 



E. 



±s^/A + B + ±t^/A-B 



s,t 



(27) 



where and ±f are two independent signs and 

A^V'^ + + t^a^ + g\c 

= {VV + tg+o?f + t^gWP^, 
a = 2cos(fc), /3 = 2sin(fc). (28) 

Taking the ground state of -ffo(3) as a trial state the variational 
ground state energy per site is (taking V — V for simplicity) 



3t 



^var,0(3) 



bonds 



bonds 



(29) 
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Here and throughout the rest of this section we will drop the 
indexes on the operators since only averages of bilinears with 
two equal indexes are non-zero because of the 0(3) symme- 
try. Minimizing this variational energy is typically equivalent 
to a mean field calculation. We have found two classes of so- 
lutions to the mean field equations that give low values of the 
variational ground state energies: one that is translationally 
invariant and one that is dimerized. 



1. Translationally invariant phase 



where a_|_ = (i + g+)/V and E{x) is the complete elliptic 
integral of the second kind.^" From this we can compute the 
translational invariant averages of the operators that appear in 
the variational calculation by differentiation with the result 



'{ill) 



dt 



E{-al) 



(31) 



Here K{x) denotes the complete elliptic integral of the first 
kind.-"' Interestingly, if we view the ground state of the mean 
field Hamiltonian as a trial state, the bound on the ground 
state energy does only depend on the variational parameters 
9+ and V in the combination a+. This means that there is 
a one-parameter family of mean field Hamiltonians that have 
the same ground state. We can imagine to try to fix the best 
value of a+ in four ways: 1) We use the mean field Hamil- 
tonian to construct a trial density matrix at finite temperature 
T. Taking the limit T ^ in the trial free energy the entropy 
term is maximized if one minimizes the gap in the mean field 
Hamiltonian. In our system this means that we should choose 
g+ — 0. 2) We can consider a calculation to second order 
in the interaction Hamiltonian Hi — Hy^usa — Ho{y)- 3) We 
maximize the energy gap in ffo(3) by taking g^ — t. This has 
the additional appealing property that every eigenvalue in (27) 
becomes double degenerate since A = B with this choice. 4) 
We use the standard Hartree-Fock decoupling procedure. 
The usual Hartree-Fock scheme gives 



g+ = 4t(i^/i)^ 



(32) 



and picks out particular values of g+ and V. Note that the 
expectation value in the atomic ground state gives V = 3 J/4 
and g+ = 0. The mean field self-consistency conditions (31) 
and (32) are easily solved numerically. The variational energy 
is exact in the limit J/t — ?► cx) but gives —t{3/TT + in- 
stead of the correct value — M/tt if the limit J/t — ^ 0. In 
the limit J = the energy of this state is therefore about 
15% too high and hence not a good approximation to the 
ground state. The trial state is better at intermediate values 
of t/J: taking J = 1 and t/J = 1/2 the best variational 
energy is e„ar,0(3) ~ —0.878, which should be compared 
to the most accurate estimate from high-order series expan- 
sions e « —0.926.^' The result is therefore about 5 % to 
high for these parameters. The discrepancy can presumably 
be made smaller by considering fluctuations around the mean 
field state. 



This phase has V = V > 0, g- =0, and g_|_ > 0. In 
this case the ground state energy per lattice site of Ha in each 
component is 



^ — (30) 



2. Dimerized phase 

For smaller values of J/t a dimerized solution, which 
has the form of a spin-Peierls state, is found to be energet- 
ically favorable. It is characterized by y = ^ > and 
g- = ±{g+ - S) with < S < g+. The state wifli (5 = 
is maximally dimerized and has all of its /-spins locked up 
into singlets with one of its nearest neighboring /-spins. As 
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5 grows the dimerization diminishes until it goes away when 
9-=0- 

For J/t ^ 1 the optimal value of g_ is found to be ex- 
tremely small and no gain in energy is found compared with 
the translationally invariant phase. For J/t <^ 1 we find that 
(5 <C .g+, but a non-zero i5 is needed for a self-consistent solu- 
tion. Taking the state with 6 — as a variational wave func- 
tion we find that this solution is energetically favorable to the 
translational invariant phase for J/t < 1.24 in ID. In the limit 
J/t ^ it gives -t{3/TT + 1/4) which is about 5% too high, 
but clearly favorable to the translationally invariant state. 



B. SO(2) X Z2 -symmetric mean field at half -filling 

In the 0(3)-symmetric mean field states there are three de- 
generate fermion bands. This leads to a natural description 
of triplet excitations in terms of two quasi-particle excita- 
tions, but it is unnatural in terms of the original description 
in terms of c-electrons and localized /-spins. The 0(3) sym- 
metry of the Hamiltonian is also broken in the presence of a 
chemical potential or crystal fields. In this subsection we will 
therefore allow for a less symmetric mean field solution that 
has a global SO(2)xZ2 symmetry with the SO(2) generator 
+ + I^. We can then write the mean field Hamiltonian 
as i?so(2) = H12 + i?3, where H3 is of the same form as in 
Eq. (23). The mean field Hamiltonian in the remaining com- 
ponents can be decomposed into an onsite part and a hopping 
part according to H12 — ^12+^12- The allowed terms in this 
Hamiltonian are restricted by symmetry. The general on-site 
local term is 

= ^(171^1 + iJ2^J■2) + V{ijifli + 172^2) 
+ mo{ifj,ifi2 + ilil2) + mo{iilijl2 + «7i72) 
+ ™i(i7iM2 + «Mi72) + n^i(«7iM2 + iAi72) 
+ m3(i/ii^2 - «7i72) + mz(ijlijl2 - «7i72), (33) 

and the general hopping term is 

HI2 = -^Y^ ^lla{ri)la{rj) 
a=ia {ij) 

+ ^Hi-Tj -ri)i^la{r^)^la{rj) 

a=l,2 (jj) 

+ X! ~ ■^»)[*Aii(n)A2(^j) - i^-2{ri)fii{rj)\. (34) 

Note that we use the same notation for some of the variational 
parameters in Hj, and H12, i.e. g and V , although their val- 
ues will in general be different. All in all there are 16 mean 
field parameters (in ID), with 4 in the third component, in the 
most general SO(2)-symmetric mean field Hamiltonian. As 
in the previous section the spectrum of H12 can be found by 
going to Fourier space using (24) and diagonalizing the result- 
ing 8x8 matrix. Generically the spectrum is then given by the 
solutions to a quartic equation for the pair (a, /3) and one for 



{—a. P), but if the quantity 

/i = (mo — mQ){m\ — toJ + 777,3 ^ ^'^s + ~ y) 

-(7773 + 7?73)(i'a2-G+-G^), (35) 

vanishes the spectrum is again given by (27) with different A 
and B. We have also introduced G+ — g+a + and 
G_ = g-p — h+a. We defer a fall analysis of the gen- 
eral mean field to a later study. In the following we con- 
fine ourselves to the simplified situation with 77ia — —rria 
for a ~ 0, 1, 3 (to allow for antiferromagnetism) and V = V 
(leading to sublattice-symmetric on-site {i'yi-i)) so that /i = 0. 
Then the spectrum is given by (27) with parameters 

A = fa^ + G\ + G'^_+ 2{V^ +ml + ml + ml), 

= {V^ + 777^ + 777^ - mlf + {mo + m^ft^a^ 

+ (Gl + Gl)[imo - 7773)2 + t^a^] 

+ 2G+{V^ ~ ml)ta + AG-miVta, (36) 

and the same with a —a. This reduces to (28) when 7770 = 
™i — ?7i3 = h± = and V = V . We also note that this 
spectrum is the same as the mean field spectrum of Ref 16 
when 7771 = h± = g± = 0. The imposed symmetries leads to 
the following relations for the operator averages 

(47172) = -(«7i72), 

(7^1^2) = -{ijllij.2), 

(*7iAii) = (472^2) = («7iAi) = («72/i2>, 
(»7iAi2) = (iMi72) = -(i7iM2) = -(iMi72), 
{i^il^j.2) = -(iM2Mi>- (37) 

The relations between the averages on different sublattices can 
be derived by considering inversion symmetry about a bond 
in the crystal together with translational invariance. The vari- 
ational energy per site is 

e™r,so(2) = - ^ X! (2(*7i7i> + («7373)) 

bonds 

bonds 

+ ^(*7iMi> + ■^(«73Ai3) - ^(«7iMi)(«73Ai3) 

- ^(*7iMi)^ - ^(«7iM2)^ + ^(j7i72)(«MiM2), (38) 

where we have used the symmetries in (37). Comparing this 
with the 0(3) case in Eq. (29) the difference is the possibil- 
ity of having nonzero averages for terms involving a coupling 
between the first and second component on the second and 
fourth line. Of particular importance is the possibility of hav- 
ing (7/11/^2) 7^ since this will allow the system to take full 
advantage of the kinetic term in the limit J 0. The op- 
erator averages can be calculated by taking the appropriate 
derivatives of the mean field ground state energy, just like in 
the 0(3) case. The expressions we need are given in the Ap- 
pendix in Eqs. (A13) and (A14). We now have all the pieces in 
place [i.e., Eqs. (27), (36), (38), (A13), and (A14)] to perform 
the variational mean field study. 
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1. Result of the S0(2)-symmetric mean field study 



to show that 



The values of mi, both 5-'s, and h± are found to be ex- 
tremely small when minimizing the variational energy for all 
values of J/ 1. Thus we will set these parameters to zero in the 
following discussion. For large values of J /t the variational 
parameters flow towards the family of 0(3)-symmetric solu- 
tions. For J /t < 1.56 nonzero 777,3 ^nd 7779 are found to lower 
the energy with respect to the 0(3) family. This implies an 
antiferromagnetic SO(2)-symmetric solution that always is fa- 
vorable to the dimerized 0(3)-symmetric solution at the mean 
field level. This is not surprising in view of earlier work in 
ID,"* 20,'"^ and 3D." In the Hmit J/t the antiferro- 
magnetic solution reproduces the correct value of the ground 
state energy, namely — 4t/7r. In our mean field analysis the 
transition to the antiferromagnetic state is discontinuous. The 
variational mean field energies for the different trial states are 
shown in Fig. 1 . 



3^ 




FIG. 1. Variational energies for different trial states in ID. From 
top to bottom: 0(3)-symmetric translationally invariant state, 0(3)- 
symmetric dimerized state, and S0(2)-symmetric antiferromagnetic 
state. Only the last state gives the correct value in the limit J/t — > 
0. The mean field theory predicts a transition between the 0(3)- 
symmetric spin liquid phase and an SO(2)-symmetric state with an- 
tiferromagnetic order at J/t ~ 1.56. 



E (i/l„(0)/Xa(r)) 

1 

'n]2 



^'sin(fc.r)(l-2(/it(fc)^Jfc))). (39) 



In the translationally invariant 0(3) phase g_ = 0. In 
this case particle-hole symmetry enforces (/ij,(fc)/i^(fc)) — 
1/2 for all fc, which means that Xtit^if) — 0. Similarly 
{ifi3{r)fj.s{r')) = in the SO(2)-symmetric antiferromag- 
netic state. When the two /i operators reside on different sub- 
lattices we have 



= ^E'c°«('^-^)<*^a(fcR(fe)>+h.c., (40) 



in the translationally invariant 0(3) phase, and from the diag- 
onalisation of (26) we obtain 



(i/xl(fc)/2jfc))+h.c. = - 



(41) 



With this result it is straightforward to evaluate the sum in 
(40) in the continuum limit numerically. The result is illus- 
trated in Fig. 2 and clearly shows that the result is an alter- 
nating almost exponentially decaying function of r. In the 
atomic hmit a+ <C 1, and the non-local correlations are small: 
Xi-ifii^) ~ ~'^+/4- This will have direct consequences for 
the spin-spin correlation functions and the triplet pairing am- 
plitudes. Similarly we can show that (77a(0)7a(r)) ~ 0, 
and {i'ya{0)ja{i')) = —X;^/: (''') in the translationally invariant 
0(3)-symmetric phase. It is interesting to note that these cor- 
relation functions depend on the variational parameters only 
through a+ and not on the actual spectrum of i/o(3)- 



V. CHARACTERIZING THE 0(3)-SYMMETRIC PHASE 

In this section we will characterize the 0(3)-symmetric 
mean field phase further by looking at the spin-spin correla- 
tion functions and the triplet pairing amplitudes. First we note 
that the average c-electron spin and charge (measured with re- 
spect to half-filling) as well as the average /-spin are zero on 
every site. As we will see, the spin-spin correlation functions 
are rotationally invariant, and we are thus dealing with a spin 
liquid state. That the system is a gapped spin liquid for large 
values of J/t is known,' but the 0(3)-symmetric trial state 
provides a simple realization of such a state for finite values 
oft/J. 

To characterize the state we first derive the correlation func- 
tions for some of the fermion bilinears. Note that we do not 
consider the trivial autocorrelation functions in the following. 
Using the Fourier representation in (24) it is straightforward 



^ -10 

I 



■15 



1 2 3 4 5 6 7 
n=(r+a)/(2a) 

FIG. 2. The logarithm of the absolute value of the correlation func- 
tion Xt^iiif) = {ilJ-aiO)p.a{i')) at the value a+ — 0.8, which is 
appropriate for the trial state at t/ J ~ 0.5. The dashed line is the fit 
to a straight line corresponding to exponential decay. 
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A. Spin-spin correlation functions 

In the 0(3)-syminetric phase the spin-spin correlation func- 
tions are easily shown to be rotationally invariant: 

{Sl,{r)Si,{r'))^S'^S^a.'{r' -r), (42) 

for all combination of a, a' = c, /. This is a consequence of 
the representations in (9), (10), (12), and the definite fermion 
parity of all three components in the 0(3)-symmetric ground 
state. It is also possible to work out explicit expressions for 
the correlation functions in detail, as an example let us look 
at / spin-spin correlation function. Using the result of the 
last subsection, Sff{r) = if the spins reside on the same 
sublattice and Sff{r) < otherwise. Explicitly we have 

{S}{0)S}{r))^-6'^^xUr). (43) 

From the behavior of X/j/i(^) we see that this is a rapidly de- 
caying negative function. 



B. Superconducting correlations 

It is easy to see that the singlet Cooper pair amplitude is 
zero in both the 0(3) and the SO(2) phases. This is a con- 
sequence of the representation in (11) and (10) and the defi- 
nite fermion parity of the third component in the mean field 
ground states. The same argument shows that the triplet pair- 
ing amplitude (c^(r)c^(r') + c^{r)c^{r')) vanishes. Using 
the SO(2) symmetry around the third axis we also find that 
{c^{r)c^{r')) — 0. The third triplet pairing amplitude in the 
translationally invariant 0(3) phase is 

A;;(r) = (c^(O)c^(r)) = 

(73(0)73(t-))/2 - 2(M3(0)M3(r-))^ (44) 

This expression vanishes in the translationally invariant 0(3)- 
symmetric phase because of the particle-hole symmetry. It 
remains the last triplet pairing amplitude on different sublat- 
tices. In the 0(3)-symmetric phase this is given by 

A^;(r) ^ (c,(0)5,(r)) 

= (*73(0)73(r))/2 + 2{tfi3iO)il3ir)f. (45) 

The behavior of A^^(r) is therefore simply related to Xma(^) 
and is a rapidly decaying alternating function. The results of 
this section clearly shows that spin-rotational symmetry, and 
hence time-reversal symmetry, is broken in a subtle way in the 
superconducting pairing correlation function in this state. 

VI. CONCLUSIONS AND OUTLOOK 

The main finding in this work is an 0(3)-syinmetric repre- 
sentation of the bipartite Kondo lattice model at half-filling. 
To the best of our knowledge this representation has not been 



written down previously. We have used this representation to 
construct and investigate an 0(3)-symmetric mean field state 
in ID, and found it to be a good trial wave functions for 
large to moderate values of J/t. At smaller values of J/t 
a state with antiferromagnetic correlations is favored. The 
0(3)-symmetric state is a gapped spin liquid with rotation- 
ally invariant spin-spin correlations and a finite (short-ranged) 
triplet pairing amplitude. 

For the future it would be interesting to apply the transfor- 
mation to other lattices than the simplest ID case considered 
here, and to see what this representation can tell us about the 
Kondo lattice away from half-filling. Another direction of re- 
search would be to study the effects of the finite triplet pairing 
amplitude and to allow for slow fluctuations in the direction 
of this order parameter. 
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Appendix A: Conventions and mathematical details 

1. Original basis 

We follow the convention of Ref. 10 to enumerate the states 
in the original basis: 

|1>-|0>, 

|2) = 4,t4Jo)' 

|5) = 4,t4,^4.t4,ii'^^' 
|8) = c|;_^4^^io), 

|9>-4,J0)' 

|io>-4Jo)' 

|ii) = 4i4.t4Jo)' 
|i2>-4t4;4,Jo)' 

|13)=4,^|0), 
|14>= 4-^10), 

|15)-4,t4.T4,J0)' 

|16)=4,t4,;4.t|0)- (A1) 



10 



2. Time-reversal symmetry 



vacuum state to \0)s- Explicitly 



We choose the phase convention for the action of the time- 
reversal operator T on a generic spin-full fermion level de- 
scribed by the operators and to be 



(A2) 



This convention implies that 1) The singlets |1) to |5) are in- 
variant under time reversal as expected. 2) The t states trans- 
forms into the corresponding I states (and vice versa) with our 
phase convention. 3) The triplets transform as 



T\l = l,m) = -(-l)"|/ = l,-m), 



(A3) 



which is different from the transformation of the conventional 
spherical harmonics.-^ Restricting ourselves to real coeffi- 
cients only the triplet (|6) + |8))/-\/2 is invariant under time- 
reversal. 



3. Unitary transformations 



Let us first define two different basis sets: is defined 



as in Eq. (Al) with c[ , 
for {IO2} with 4^ c, 



-I 



c2,a and c]^^ c}2,^. 



Similarly 
The uni- 



tary operator U implements the transformation between the 
old basis {|Z)i} and the new basis {IO2} via the relations 



\l)2 = U\l), 



for / = 1 ... 16. 



(A4) 



In the following we will define U to be the matrix representing 
U in the old basis, i.e. 



Uk,i = i{k\U\l)i = i{k\l)2. 



(A5) 



From the definition of the states, and completeness of the ba- 
sis, we see that the creation operators in the two bases are 
related by 



4 



(A6) 



This clearly preserves the fermionic anticommuation rela- 
tions. We also have c/2,(t|0)2 — iff c/i,cr|0)i ~ etc., 
which is consistent with the notion that annihilation operators 
annihilates the vacuum. 



4. The electron-hole basis 



The electron-hole basis of Ref. 10 is obtained by substitut- 
ing h\ and cj, ^ — ej^ in Eq. (Al), and changing the 



1) eh 

2) eft 



= |0)., 



e|et|0). 



3) eh 

4) ,„ = h\hl\Q),, 

5) eh 



{h\el-h{e\m 



I eh — 



h\hle\el\Q), 
1 



7) eh = -j^{h\el + 

8}eh = h\e\\0)s, 
9)eh = e^|0)., 

|lG)eh = /l||0>s, 

|ll)eh = hle\el\0)s, 
\12),h^h\hlel\0),, 
|13)eh = e||0)„ 
\U),h^h\\0),, 

1 15) eh = /i|e|e||0)s, 
\16),h = h\hle\\0),. 



hle\)\0)s 



(A7) 



5. The electron-hole transformation 

Diagonalizing Eq. (2) and using the procedure of Sec. 11 A 
to assign the states we obtain the matrix Ueh that implements 
the transformation from the original basis in (Al) to the new 
one in (A7). Exphcitly 



Ueh 



^Us 

I3 

Ui 

C/^y 



(A8) 



with submatrices 



Us = 







/ V. _ 

sin((^2)/\/2 1/72 

COs((y92) U 

sm{(p2)/V2 -l/%/2 
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_\ 

-cos(.^2)/\/2 
sin(.^2) 
)s(^2)/V2 



- cos 







u^ = u^ 





cos(</?i) 
sin((/3i) 

sin((pi) cos((^i 

\ cos((/7i) sm{(p 




J 
(A9) 



(AlO) 



I3 is the 3x3 unit matrix, tan((^i) 

+ 4H/2)^ and tan(^2) = 41¥/(C/ + ^/TPTwm). 



2W/{U 
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6. Basis for the low-energy sector in the Kondo limit 

To describe the Hilbeit space in the low-energy sector in 
the limit [/ — > oo we only need three operators. It is therefore 
natural to use one operator (i.e. c\) to describe excitations in 
the high-energy sector. We use the following convention to 
label the states 

Wk = \0)s, 
\2)k^c\\0},, 

|3)K = 4|0)a, 

|4)i^ = 4|0)„ 

|5) if = 44410),, (All) 

|6)if = 44l0)s, 
|7)^ = 44|0)„ 

|8>/f- 44 |o)s, 

\m)K = 4(-l)"'^"'^"' \m - 8) if, m = 9 ... 16. 

If we are only interested in the low-energy sector of the theory 
the assignment of the states in the high-energy sector does 
not matter. We can therefore make any convenient consistent 
choice, for example 

\9)K = \5)eh, |10)K=|ll>eh, |11)a' = |12)e/., 
|12)a' = |15)e/., |13)K = |16)eh, \U)K = \2)eh, 
|15)a- = |3)eh, \l6)K = \'^)eh- (A12) 



If the high-energy operator c\ is of interest one should make 
a better informed choice. 



7. Expressions for the operator averages 

The expressions that we need to calculate the operator av- 
erages in the SO(2)-symmetric mean field states are 

X] (*7373> = -29(63, 
bonds 

(^73^3) = dve3, 
(W3(±l)) = (5g+±93Je3, (A13) 

for averages in the third component and 

^{■ili'li) = -<9tei2, 

bonds 

2(z^i/ii(±l)) = {dg^ ±dg_)ei2, 

2(«7iMi> = 5yei2, 

2(i^l/i2) = {d,no + 5,„3)ei2, 

2(j7i72) = {dmo - 5m3)ei2, 

2(^71^2) = 9„,ei2, (A14) 

for averages in the other two. 
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